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Abstract. The spin one-half Heisenberg chain with U q [SU(2)] symmetry is studied via density-matrix 
renormalization. Ground-state energy and g-symmetric correlation functions are calculated for the non- 
hermitian case q — exp(i7r/(r + l)) with integer r. This gives bulk and surface exponents for (para)fermionic 
correlations in the related Ising and Potts models. The case of real q corresponding to a diffusion problem 
is treated analytically. 



1 Introduction 

The density-matrix renormalization method (DMRG) in- 
troduced by White |lj has opened new possibilities in the 
study of one-dimensional quantum systems. With this ap- 
proach not only the energies of large systems can be cal- 
culated very accurately, the method also yields correla- 
tion functions. M. It was tested e.g. on Heisenberg |Q and 
transverse Ising j3j chains and has been applied to a vari- 
ety of systems like higher-spin chains Q , ladders || , Hub- 
bard chains impurity models [Q, phonon systems Q 
or transfer matrices (IJlij . It is therefore natural to try it 
in still other situations and to see how it performs there. 

The system which we have studied is the U a \SU(2)]- 
symmetric Heisenberg chain with Hamiltonian |TJ] 
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with open boundary conditions and 
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where the erf are Pauli matrices. Some basic features of 
the g-symmetry are sketched briefly in an appendix. 

The motivation for considering this problem was two- 
fold. Firstly, the model comprises various other physical 
systems and situations which thereby can be treated in a 
unified way. In particular, we were interested in calculating 
correlation functions and their exponents. Secondly, the 
Hamiltonian is non-hermitian for complex q and offers a 
possibility to investigate this situation. 
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In the antiferromagnetic case (sign — ) and for complex 
q = cxp(i7r/(r + 1)), r integer, the chain is related to the 
critical models with central charge c = 1 — 6/(r(r + 1)), 
in particular the Ising (r = 3) and the three-state Potts 
model (r = 5) [|l2|]l3[l . For certain real values of q it corre- 
sponds to Potts models with more than three states p4[ . 
The ferromagnetic chain with real q, on the other hand, 
describes a non-equilibrium problem, namely the hopping 
of classical hard-core particles on discrete sites, with a bias 
in one direction [p^[ . For q = 1, finally, one recovers the 
simple isotropic Heisenberg model. 

The Hamiltonian (Q) has real eigenvalues and can be 
treated by Bethe ansatz. In this way, finite-size spectra 
have been obtained in connection with conformal invari- 
ance . Using the Temperley-Lieb properties of the op- 
erators e,i, one can also find a representation by hermitian 
IRF operators, for which some DMRG calculations have 
already been done []l6|| . However, as mentioned above, we 
are interested in non-hermitian effects. These can be var- 
ied in this model by changing the parameter q. 

The general arguments for the DMRG procedure Q| 
do not involve the Hamiltonian and the approach is there- 
fore applicable here. In fact, non-hermitian operators have 
been treated occasionally 1 10 17],|l8| . There is some free- 
dom, however, which density matrix to choose then. In 
our case we always used a hermitian one, constructed from 
the right eigenvectors of ([!]) . One then finds that its spec- 
trum shows a marked and systematic variation with r. 
The eigenvalues decrease more slowly than for the simple 
Heisenberg chain. Nevertheless, there is no difficulty in 
obtaining the ground-state energy of H to high accuracy. 
This is discussed in section |[ 

Our main effort involved correlations. The functions 
which we studied were cy-symmetric generalizations of the 
quantity < ener m > and correspond to fermionic and para- 
fermionic correlators in the Ising and Potts case, respec- 
tively ]l9| . In a previous study with exact diagonalizations 
of short systems no conclusive results could be obtained 
for the Potts model [E(|. From our calculations up to 100 
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sites we were able to obtain functions with clear scaling 
behaviour from which bulk and surface exponents could 
be extracted. The results are presented in section ||. In 
the Ising case, where one has the exact solution as a test, 
they are quite good. In the Potts case, there are still some 
uncertainties, since the bulk exponents of the four func- 
tions which we considered differ while their spinor char- 
acter suggests that they are equal. 

The diffusion problem, treated in section |], has a very 
different but also interesting nature. Here one finds that 
the eigenvalues of the density matrix drop to exactly zero 
very rapidly, so that the DMRG procedure terminates af- 
ter a certain number of states. This can be understood 
from the nature of the ground state which is known here. 
From it, the density matrix can be derived and thus one 
has a non-trivial example of a system where the DMRG 
procedure can be carried out analytically. 

A summary and some additional remarks can be found 
in the concluding section ||. 

2 DMRG procedure and ground-state 
energies 

In the DMRG approach, a selection of relevant states is 
made, not for isolated, but for already interacting blocks 
which form the total system. In our calculations we always 
used the infinite-size variant of the algorithm H . 

As usual the chain was divided in two parts labeled by 
the indices 1 and 2, respectively. Then the right ground- 
state vector \<P r >= X^ijl* >i \j >2: where \i >i and 
\j >2 denote the basis states in the two parts, and the cor- 
responding eigenvalue E were calculated by using a com- 
bination of a vector-iteration procedure (power method) 
and a Lanczos procedure with modified scalar product. 
The whole calculation was done in the S z = subspace 
which reduced the numerical effort significantly. 

From the ground-state vector the hermitian density 
matrix 

p=|* r X* r | (3) 
of the chain was constructed where 

<<P r \ = |# P > f (4) 

is the usual hermitian conjugate. From (||) the reduced 
density matrix p\ of part 1 of the chain follows as 

Pi = J2®ij$i>j\i>n<i'\ (5) 

and similarly for p 2 . Since the boundary terms in (Q) break 
reflection symmetry one has to work with both matrices 
here. 

The use of hermitian density matrices is already sug- 
gested by the considerations in jjj where p\ enters via 
an optimization procedure for a vector like \<P r >. Non- 
symmetric p's have been used in some cases P, |io|]i"g| ] but 
lead to obvious problems if their eigenvalues become com- 
plex. We tested the choice p = \<P r >< 4>i\, where < #j| is 



the left eigenvector of H, on the XY-analogue of ([jj) but 
it gave unsatisfactory results, in particular large errors 
which did not decrease properly with m. 

The m eigenvectors |A„ >i, |A„ >2 of p\ and p 2 cor- 
responding to the largest eigenvalues A„ are used as new 
basis states for the two parts of the chain in terms of which 
all necessary operators are expressed. After inserting two 
sites between part 1 and part 2 the procedure is restarted. 

For the calculation of expectation values 

<6>=<$i\6\$ r > (6) 

one also needs the left eigenvector. In the present case 
this is just the transpose of |<5 r > since H is complex 
symmetric. 

For the efficiency of the DMRG procedure it is impor- 
tant that one can calculate accurate results with a rela- 
tively small number of states of the Hilbert space. There- 
fore the eigenvalues of the reduced density matrix should 
decay exponentially. However, in general it is not a priori 
clear if the spectrum exhibits such a behaviour and one 
has to check it numerically. Figure |l| shows such spectra 
for different values of the parameter r. One can see roughly 



logX a 










-2 




r=3 

r=5 " 


-4 




r=°° 


-6 
-8 






-10 






-12 






-14 






-16 







10 20 30 40 50 60 70 80 90 100 



n 

Fig. 1. Eigenvalue spectrum of the reduced density matrix pi 
for different values of the parameter r and chain length L = 30, 
calculated with 64 kept states. 

exponential behaviour with a fast initial and a slower fi- 
nal decay. However, the results for small r lie considerably 
above those for the isotropic Heisenberg model. The q- 
symmetric chain is therefore numerically a less favourable 
case. This feature is essentially due to the imaginary 
boundary terms. When added to an XX-chain they shift 
the density-matrix spectrum in the same way. 

We have calculated the ground state energies for dif- 
ferent values of the parameter r and for different numbers 
of kept states m. Table [l] shows some results for m = 128. 
(We note that the results given here are calculated with- 
out the constant — (L — l){q + q~ v ) / 4 in ([!]).) The accuracy 
can be checked directly for r = 3 by comparing with the 
analytic result |2l]] . In this case the deviations are less than 
1 -9 which is in accord with a truncation error of the or- 
der of 10 -9 in the density matrix. One can also compare 
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L 


r=3 


r=5 


4 


-0.6604497504 


-0.7393376252 


10 


-0.7443482698 


-0.8012448364 


20 


-0.7736274980 


-0.8237779981 


40 


-0.7885618957 


-0.8354921842 


60 


-0.7935870982 


-0.8394707437 


80 


-0.7961088422 


-0.8414746652 


L 


r=7 


r=oo 


4 


-0.7687664240 


-0.8080127018 


10 


-0.8227289432 


-0.8516070414 


20 


-0.8427349544 


-0.8682473334 


40 


-0.8532195604 


-0.8770736642 


60 


-0.8567946053 


-0.8801004992 


80 


-0.8585981414 


-0.8816309100 



Table 1. Ground state energies per spin for different r, cal- 
culated with m = 128 for r = 3, 5, 7 and m = 64 for r — oo. 



with the values obtained by Sierra et al. [[L6| in a different 
DMRG calculation with rn = 160. We could reproduce 
their data up to 9 digits. 
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Fig. 2. Ground state energies per site for different values of 
r calculated with m = 128 (for r — oo with m = 64). The 
symbols denote data from Sierra et al. |16[. 



Plotting the ground-state energies per spin vs. 1/L, 
one obtains the well-known asymptotic behaviour 



E 

T 



B 

T 



c_ 

L 2 



(7) 



where B is the boundary energy and C measures the 
Casimir effect. This is shown in figure || for several val- 
ues of the parameter r. 

One expects that for a non-hcrmitian operator the pro- 
cedure does not necessarily give an upper bond for E. This 
is in fact the case. For r = 3, the deviation from the exact 
result, although quite small, changes sign for a certain L. 
Similarly, the approach of E towards its finite-size limit 
as a function of m is not monotonous. This is shown in 
figure ||. One notes, however, that the behaviour becomes 
monotonous for large r. 



E/L 




-0.838 



Fig. 3. Ground-state energies per site for r — 3 and r — 5 vs. 
the number of kept states m of the density matrix for chain 
length L — 50. 



3 Correlation functions 

For the (/-symmetric chain with q a root of unity, the usual 
spin correlation functions are not the appropriate quan- 
tities. Instead, we consider operators which are adapted 
to the U q [SU(2)] symmetry. In this way one is lead to 
generalizations of the scalar product cricr rn of two spin 
operators given by p9| 



dtrr 



1 



-l,m 



(8) 



with 



,±1 



(9) 



where the ej are defined in equation (^). Written out ex- 
plicitely, these are non-local objects containing strings of 
spin operators. However, they are constructed solely in 
terms of the ej from which the Hamiltonian (Q) is built. 
Therefore they can be translated into any other represen- 
tation of the Temperley-Lieb algebra which the ej fulfill 
[ p^[ . In the cases r = 3 and 5 such representations are 
given by the critical transverse Ising model and the Potts 
model with L/2 sites and open boundary conditions. In 
these cases < gf m >~< gf m >, thus we work in the fol- 
lowing only with one operator and drop the superscript. 

In the Ising case one finds that the string operators 
cancel the Jordan- Wigner factors in a fermion picture and 
the < gi, m > become simple fermionic two-point functions 
[ ^9| . For example 

<g 2j -i,2 k > = --^=(-iy +k 



< 92j,2k-l > 



1 

'V5 



(4+d,)(4-d fc )) (io) 

(-iy +k 



x((4-4,) (4+^)) (ii) 
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where the dj, are fermion operators. These functions 
can be calculated exactly and are given as sums for fi- 
nite L. The two quantities < <?2j,2fc > and < .92^-1. 2fe-i > 
vanish identically. The other two correspond to different 
behaviour at the boundaries where ( p~0| ) remains finite and 
( |TT| ) vanishes for fixed k and finite chain lengths. Due 
to the open ends, one is dealing with surface critical be- 
haviour here. After performing the continuum limit one 
can determine the critical exponents, and one finds the 
bulk exponent x = 1/2 for both correlation functions and 
the surface exponents x s — \/2 and x s = 3/2 for ( fL0| ) and 
(|n[), respectively |§. 

For the Ising case one can also construct the conven- 
tional spin correlation function < erf > from the quan- 
tities ej appearing in H, using the property (<t x ) 2 = 1. 
However, this is not possible for the Potts model. Also or- 
der parameter profiles, as determined by a direct DMRG 
calculation in P2f , cannot be obtained in the present for- 
mulation. However, the g's are also interesting quanti- 
ties since they become parafermion operators in this case 
]l9| , p3] | . Their exponents will be studied in the following. 

We considered the correlation function < gix/2 > 
(i = 1, ...,L/2 — 1) with one point fixed in the mid- 
dle of the chain while varying the position of the other 
one. In this way profiles were computed for chains with 
up to 100 sites. The correlations were significantly less 
accurate than the energy, so that in general m — 200 
states were kept. This was a practical limit with respect 
to storage and time. A complete run for fixed r took sev- 
eral weeks CPU time on middle-performance DEC- Alpha 
workstations. The amount of memory was about 250-300 
MB RAM and up to 500 MB harddisk space. 

An example of the resulting correlation function is 
shown in figure |J. The oscillations are due to the anti- 
ferromagnetic character of the chain. Except for the case 
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Fig. 4. Correlation function < <?;.l/2 > f° r r = 5 and chain 
length L = 82, calculated with 200 kept states. 

r = 3, where only the two functions ( |l0| ) and ( pd| ) are non- 
zero, one can distinguish four different functions given by 
the maxima and minima for even and odd L/2. 

As a check on the accuracy, one can compare with 
exact results for L = 24 ^4|. One then finds deviations 



less than 10~ 4 in the Ising case and less than 10 -5 in the 
Potts case while the truncation error is only 10~ 14 — 10~ 15 . 
The correlations are therefore much less precise than the 
ground-state energy. One also finds that the values can be 
above or below the exact result. The errors increase, as 
usual in the DMRG, as one moves towards the boundary. 
This effect can also be seen in figure S, where two correla- 




0.0001 1 1 ■ — 1 ■ — 

1 10 
i 

Fig. 5. Log- log plot of the correlation function < > for 

even i and chain length L — 82 (r — 5) and L — 90 (r = 3), 
calculated with 200 kept states. The curve denotes analytic 
results for r = 3. 

tion functions are shown for large systems. Up to R ~ 30, 
where R = L/2 — i is the distance from the center, the 
numerical values for r = 3 reproduce the analytical ones 
quite well. Near the boundary, where the function is very 
small, however, there are considerable deviations and the 
expected power-law behaviour g ~ R Xs ~ x , with R s = i 
denoting the distance from the surface, breaks down. This 
is different for r = 5, where the curve remains linear also 
in this region. One therefore expects that the results for 
larger r are more accurate. 

In the middle of the chain, the situation is more favour- 
able since the quantities have undergone fewer iterations 
here. For small R (R <C L/2) one finds the expected bulk 
behaviour g ~ R~~ 2x , as shown in the log- log plot of figure 
^. However, if one determines the exponent x from this 
data for r = 3, one still has a relative error of about 5%. 

In order to improve the accuracy, we used the scaling 
form for g 

where F(R/L) — > const for R/L — > 0. Collecting results 
for different L we constructed a scaling plot R 2x g vs. R/L 
and varied x until the data fell onto a curve as well as 
possible. The result of such a procedure for r = 5 is shown 
in figure |?]. It proves that scaling is indeed fulfilled, in 
contrast to |2(J where it could not be seen for chains up 
to L = 28. The exponent x for r = 3 then is only about 
2% off the exact result x = 1/2. 

A further improvement could be achieved by extrapo- 
lating the bulk values of g to the infinite-chain limit. Then 
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Fig. 6. Log-log plot of the correlation function < gn t L/2 >, 
R = L/2 — i, for odd R and chain length L = 82 (r = 5) and 
L — 90 (r = 3), calculated with 200 kept states. The curve 
denotes analytic results for r — 3. 
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Fig. 7. Finite-size scaling behaviour of the correlation function 
< <7h,l/2 >, R — L/2 — i, for r = 5 and five chain lengths from 
L — 66 to L — 82, calculated with 200 kept states. 



a linear fit of the data in the corresponding plot gave bet- 
ter Ising exponents. In this way, the results in table ^ were 
determined. The confidence interval of the fit parameters 
would give rise to an absolute error in the exponents less 
than 10 -3 , but the error in the data points is not known. 
Thus we cannot give a reasonable error bar but we can 
use the results for r = 3 as a clue. 



r 


< 92j,2k > 


< <?2j,2fc-l > 


< <?2j-l,2fc > 


< <?2j-l,2fc-l > 


3 




0.499 


0.502 




4 


0.386 


0.521 


0.513 


0.406 


5 


0.397 


0.520 


0.508 


0.409 


6 


0.392 


0.516 


0.502 


0.417 


7 


0.412 


0.499 


0.484 


0.422 



Table 2. Critical exponents x of the four correlation functions, 
calculated with m = 200. Exact results are known for r = 3 
and r — oo where x = 1/2. 



The final Ising exponents differ only by 10 -3 from the 
exact value and thus are very good. The situation is more 
complicated, though, for other cases. The Kac table [^5| 
for r = 5 contains 10 different conformal dimensions and 
allows for a large number of combinations x — A + A. 
However, if one fixes the spin s of the operators to simple 
values, only x = 11/20 (s = 1/2), x = 7/15 (s = 1/3) and 
x = 2/3 (s = 2/3) remain. These values were also obtained 
directly by relating the three-state Potts model to a Gaus- 
sian model p6[ . From the explicit form of the parafermion 
operators (products of local and string operators) [[l9] , it 
can be seen that they have spinor properties as discussed 
in |p6| with s = 1/3. Thus one would expect a common 
exponent x = 0.467 for all four g's, whereas the numerical 
values lie about 10% higher and lower. In principle, they 
would allow also other a;-values connected with strange 
spin. A similar situation is found for r = 4 (tricritical 
Ising model) where x = 7/16 (0.438), x = 19/40 (0.475) 
and x = 43/80 (0.537) are close to the measured values. 

The reason for the discrepancies is not clear. As noted 
above, one would expect better accuracy for the Potts 
model. If the re-values actually differ, it could perhaps be 
due to a cancellation of leading singularities in the g's 
which, when written in the Potts representation, are sums 
of two contributions. The other possibility is that in spite 
of the satisfactory scaling behaviour the numerical results 
are not yet good enough. 

We also determined the surface exponents x s . This was 
done with the help of another scaling plot. For i — R s , the 
scaling form of g can be written as 



R! x g(R s ,L) = 



X s -\-X 



G 



(13) 



where G — * const for R s /L — > 0. Using the bulk expo- 
nents as determined above, one can now find the surface 
exponent by tuning x s until the scaling behaviour ( |l~3] ) 
is fulfilled. In this procedure the first few sites near the 
boundary, for which the errors are large, were left out. 

For the Ising case we found x s — 1.58 for the func- 
tion ( fio| ) and x s — 0.56 for the function (p"i"|). One notes 
that these results differ more from the exact values than 
the bulk results. For the Potts case the following expo- 
nents were found: x s — 0.81 for < <?2j,2fc >, x s — 1.05 for 
< ff2j,2fc-i >, %s = 0.63 for < Q2j-i,2k > and x s = 0.55 
for < g 2j -i,2k-i >■ 

We have also determined the exponents for higher r. 
The values tend to their Heisenberg limits x = 1/2 and 
x s = 1 for r — > oo. However, while the exponents of the 
function < 02j,2fc— l > approach this limit from above, 
the exponents of the other functions always lie below it. 
Furthermore, the limiting values are not easy to obtain 
due to the logarithmic corrections which occur for q = 1 



4 g-symmetric driven diffusion 

It is well-known that the isotropic ferromagnetic Heisen- 
berg model also describes (in all dimensions) the hopping 
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of classical particles on a lattice with the exclusion of dou- 
ble occupancy 1 28 . This follows by expressing the master 
equation for the probability vector \P > describing the 
system 

^-\P>=-H\P> (14) 

in a spin one-half language where a spin up (down) corre- 
sponds to an occupied (empty) site. The stationary state 
is therefore the ferromagnetic ground state in the corre- 
sponding sector of fixed S z . 

A similar result holds if there is a bias in the hopping. 
In one dimension and after a canonical transformation the 
time-evolution operator then takes exactly the form (Q) 
with q given by q = (ot-/a+) 1 / 2 where a+ (a-) is the 
hopping rate to the right (left) [|l5|. The stationary prop- 
erties of this diffusion problem have already been studied 
p9f . Due to the bias, the particles accumulate at one end, 
thus producing a non-trivial density profile. In the mag- 
netic language, one is dealing with a uniaxial ferromagnet 
to which opposite boundary fields are applied at the two 
ends. These are real and H is hermitian in this case. 

The numerical treatment leads to a density-matrix 
spectrum as shown in figure [| The eigenvalues decrease 
faster than exponentially and become effectively zero be- 
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Fig. 8. Eigenvalue spectrum of the reduced density matrix pi 
of the g-symmetric hopping model at half-filling for different 
values of the parameter q, calculated with 16 kept states and 
chain length L = 30. 



yond a certain point depending on the size of the system. 
The decrease also depends on the asymmetry parameter 

q- 

To understand this, it is useful to consider the case 
q = 1 first. Then the ground state of the chain is a spin 
state \S,M > where S = L/2 and M = S z can be chosen. 
Taking M = and dividing the chain into two halves, it 
can be written as 

+s 

\<P >= \S,Q >= c m \s,m> 1 \s,-m> 2 (15) 



where \s,m >i,2 are the corresponding spin states of the 
two subsystems with s = S/2 = L/A and the c m are the 



appropriate Clebsch-Gordan coefficients. Thus \<P > is the 
superposition of (S + 1) terms and the density matrix 
pi of part 1 of the chain has only that many non-zero 
eigenvalues. It reads explicitely 



Pi 



E 



\s, m >i c m !< s, m\ 



(16) 



so that these eigenvalues are given by (? m . 

For q ^ 1 the situation does not change, because there 
are analogous decomposition formulae for the g-deformed 
angular- momentum states Bfl . In this case one finds 



1 



[2,s] 



I 2 



i/ps]! [s - to]! [s + m] 



with the notation 



T 



(17) 



(18) 



and 



[n]! = [l]-[2]....>]. (19) 
reduces to n and (|l^) becomes the normal 



For q — > 1, [n 

Clebsch-Gordan coefficient. One can check that the c" m 
given by (|l7]) are precisely the density-matrix eigenvalues 
obtained in the numerical procedure. 

Due to the g-factor in ( |l7j ) there is no symmetry be- 
tween to and —to for q ^ 1 and the degeneracy of the 
eigenvalues is lifted as seen in figure ^. Specifically for 
large q one has 

Cm ~ q-( s -™) 2 (20) 

so that the coefficient c s ~ 1 dominates and \<P > becomes 
approximately a product state 



Is, S >i \s, -s > 2 



(21) 



Physically this means that one half of the chain is almost 
filled with particles and the other one is almost empty, as 
expected for large bias. Expressions for the density profile 
were given in p9| and are easily reproduced by the DMRG 
calculations. Similarly, the treatment can be extended to 
M 7^ 0, i.e. to systems which are not half-filled. Then the 
number of terms in \<P > is even smaller, namely (S + 
1 — \M\), and becomes one for a completely full or empty 
system. The DMRG procedure gives exact results once m 
exceeds this value which increases at most linearly with 
the system size. 



5 Conclusion 

We have studied the g-symmetric Heisenberg chain for 
various cases of physical interest. For the case of com- 
plex q we showed that the DMRG procedure works well 
even though the density-matrix spectra are less favourable 
due to the non-hermitian boundary terms. We found the 
ground-state energies to high accuracy and calculated also 
generalized correlation functions, for which we determined 
the critical exponents. 
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The study of various different models via the spin one- 
half operator H was also motivated numerically, since the 
dimension of the matrices in the DMRG procedure is lower 
then. The results show, however, that to some extent this 
advantage is compensated by the sensivity of the correla- 
tions in this formulation. One should also note that the 
equivalent Potts chain has only half the length of the q- 
symmetric chain. Thus it seems that one cannot really cir- 
cumvent certain features of the corresponding problems. 
To improve the results further, one would have to increase 
the numerical effort, by using larger values of m or the 
sweeping procedures in the finite-size algorithm jlj. 

The ferromagnetic chain for real q was seen to have 
very different features and is interesting in other respects. 
Firstly, as an example where the density matrix can be 
found analytically and the DMRG procedure automati- 
cally gives the exact result. There are only a few other 
cases of this kind, for example two coupled oscillators jn| 
or finite-dimensional matrix-product states p2| . Secondly, 
it describes not only a magnetic problem, but also a non- 
equilibrium system. The use of the DMRG in this field, 
where the time evolution operators are in general non- 
hermitian, has only started [l8| and further interesting 
applications can be expected. 



Appendix 

The q- or quantum group symmetry (see e.g. fl3|| ) is a 
generalization of rotational symmetry and characterized 
by a modified commutator 



[S + ,S~] = [2S Z ] 



(A.1) 



between the generators. Here the bracket is defined in ([18]) 
so that a g-dependent function of S z appears on the right. 
For q = 1, the usual angular momentum algebra is recov- 
ered. For a chain of spins, S z — E°n/2 nas the usual 
form, but S' ± are given by 



1 L 

n=l 



E 



4« 



E <i 2 



(A.2) 



The property [H, S f± ] = of the Hamiltonian ([!]) can 
be used, as in the case q = 1, to obtain all ferromagnetic 
ground states from the one with all spins up via repeated 
application of S~ (29). In this way one can also derive the 
Clebsch-Gordan coefficients (llT 
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